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Quantum electrodynamics in 2+1 dimensions is an effective gauge theory for the so called algebraic 
quantum liquids. A new type of such a liquid, the algebraic charge liquid, has been proposed recently 
in the context of deconfined quantum critical points [R. K. Kaul et al., Nature Physics 4, 28 (2008)]. 
In this context, we show by using the renormalization group in d — 4 — e spacetime dimensions, 
that a deconfined quantum critical point occurs in a SU(2) system provided the number of Dirac 
fermion species Nf > 4. The calculations are done in a representation where the Dirac fermions are 
given by four-component spinors. The critical exponents are calculated for several values of Nf. In 
particular, for Nf = 4 and e = 1 [d = 2 + 1) the anomalous dimension of the Neel field is given 
by 7/jv = 1/3, with a correlation length exponent v = 1/2. These values change considerably for 
Nf > 4. For instance, for Nf — 6 we find r)N ~ 0.75191 and v ~ 0.66009. We also investigate the 
effect of chiral symmetry breaking and analyze the scaling behavior of the chiral holon susceptibility, 
G x (x) = (^(ar)V»(ar)^(0)^(0)). 



PACS numbers: ll.10.Kk, 71.10.Hf, 11.15.Ha 
I. INTRODUCTION 

In very simple terms we can define spinons as ele- 
mentary excitations carrying spin but no charge. They 
emerge in two-dimensional Mott insulators as a conse- 
quence of electron fractionalization. 1 Similarly to quarks 
in QCD, spinons do not emerge so easily out of a confined 
state. As a consequence, explicit observation of spinons is 
very difficult. Quarks, though being confined, are essen- 
tial to explain the properties and structure of matter. In 
a similar way, spinons, in spite of their confinement, seem 
to be essential in explaining all the remarkable properties 
of Mott insulators, including the doped ones, opening the 
way for the understanding of high-T c superconductors^ 

In several systems spinons may deconfine at a quan- 
tum critical point X One prominent example is the criti- 
cal point separating different types of order in quantum 
antiferromagnets, like for example in the phase transi- 
tion between a Neel and a valence-bond solid (VBS) 
stated which is a paramagnetic state breaking lattice 
symmetries. Spinon deconfinement also occurs in some 
paramagnetic Mott insulating states where no symme- 
tries are broken, like for example the spin liquid state 
in pyrochlore antiferromagnets;^ and in some theories 
for the underdoped state of high-T c superconductors.- 
A remarkable property of deconfined spinons is that 
they are interacting at low energies (i.e., they are aymp- 
totic interacting in the infrared), in contrast with decon- 
fined quarks, which are free at high energies (asymptotic 
freedom) £ Therefore, the field theory of a deconfined 
quantum critical point is a highly nontrivial conformal 
field theory (CFT)£ This deconfined quantum critical- 
ity (DQC) arises due to a destructive interference be- 
tween the Berry phases of a quantum antifcrromagnet 
and the instantons. This mechanism makes instanton 
events irrelevant at large distances, which in turn allows 
the spinons to deconfine. Furthermore, since deconfined 
spinons are strongly interacting, the critical exponents 



are nontrivial. Interestingly, the critical exponents fol- 
lowing from such a theory have unusual values with re- 
spect to those arising from the Landau-Ginzburg- Wilson 
(LGW) paradigm of phase transitions^ The difference 
can already be seen in mean-field theory. It is easy to 
see that mean-field theory based on the LGW paradigm 
gives an anomalous dimension ryjv = for the Neel field 
Si — (— lySxii, where i is the lattice site, S is the spin, 
and nf = 1. However, this is not the case for a mean- 
field theory in the DQC scenario. There the direction 
field is written in terms of the spinon fields using a CP 1 
representation, i.e., n^ a = ztcr a z,*, a — 1,2,3, where 
Zj = {zn,Zi2), and <r a are the Pauli matrices, and the 
constraint nf = 1 implies zjzj = 1. In mean- field the- 
ory the dimension of fields is simply given by dimensional 
analysis of the action functional. In such a case this gives 
simply dim[zi Q ] = (d — 2)/2, where d is the dimension 
of the spacetime. Therefore, mean-field theory leads to 
dim[nj, a ] = (d — 2 + t)n)/2 = 2dim[z Q ,] ) which implies 
r/N = d — 2, or r)N = 1 in 2 + 1 dimensions. Thus, we 
see that the field theory of deconfined spinons leads to a 
large anomalous dimension of the Neel field already at the 
mean-field level. Quantum fluctuations typically reduce 
this value, but it still remains considerably larger than 
the prediction of the LGW scenario**^ whose fluctuation 
corrections lead to an anomalous dimension t]n 0.03. 11 

A deconfined quantum critical point governs a second- 
order phase transition between different quantum states 
of matter at zero temperature, the paradigmatic situa- 
tion in this case being the transition from a Neel state 
to a VBS. Thus, it is important to check whether the 
proposed models for DQC actually undergo a second- 
order phase transition with a large anomalous dimen- 
sion. There has been in recent years a lot of activity in 
this direction. Indeed, some of the proposed models have 
been checked extensively in large scale Monte Carlo (MC) 
simulations 3 i 10 ' 12 i 13 ' 14 i 15 ' 16 Particularly important are 
the MC simulations for the well studied S — 1/2 quantum 
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antiferromagnet with easy-plane anisotropy. Due to its 
global U(l) symetry and the additional local U(l) sym- 
metry, the model is self-dual.— This model was predicted 
to exhibit a deconfined quantum critical point. 5,9 How- 
ever, recent MC simulation a 12 ' 13 have clearly shown that 
in this case actually a first-order phase transition takes 
place. This result is confirmed by a recent renormaliza- 
tion group (RG) analysis.- 7 Remarkably, despite the ab- 
sence of quantum criticality, the destructive interference 
mechanism between the Berry phases and instantons was 
shown to workfi^ so that we can still talk about spinon 
dcconfincment. The situation here is reminiscent of the 
Z2-Higgs theory in three spacetime dimensions, though 
in this case we have a discrete gauge symmetry rather 
than a continuous one. This model is also self-dual, 18 
but is known to undergo a first-order phase transition 
along the self-dual line.— However, the Z2-Higgs theory 
has the peculiarity of having a deconfined phase bounded 
by two lines of second-order phase transitions that be- 
come first-order when they are very near each other and 
to the self-dual line. It is not known at present whether 
a similar behavior for the easy-plane frustrated quantum 
antiferromagnet is possible. 

For the SU(2) case, a MC simulation on a Heisenberg 
model with hedgehog suppression 9 seems to support the 
DQC scenario. Another interesting model is the S = 1 
Heisenberg antiferromagnet with a biquadratic interac- 
tion between the spins ^ where numerical evidence for 
DQC has also been found. 14 Strong evidence for DQC has 
been reported in Ref. [ToL where a model featuring a four- 
spin interaction was simulated. There the anomalous di- 
mension of the Neel field was found to be t]m « 0.26. 
Further MC studies^ confirm the analysis of Ref. [lCj, in 
which a transition from a Neel state to a VBS occurs. 
There the obtained value of the anomalous dimension is 
t]n ~ 0.35. However, there is a recent paper— on the 
same model where MC simulations are reported to lead 
to a first-order phase transition. 

Due to the inherently non-perturbative character of de- 
confined quantum critical points, purely analytical and 
well controlled studies are not easy to perform. For in- 
stance, the perturbative RG applied to the SU(N) case 
can only access a quantum critical point if TV > 182.9. 17 
The smallest value of N producing a critical point, N — 
183, leads to an anomalous dimension r/N = 609/671 ss 
0.9076. In principle the result of this RG analysis would 
mean that the phase transition for the SU(2) case is a 
first-order one. However, it was argued in Ref. [Tt] that, 
similarly to the Ginzburg-Landau (GL) superconductor 
with N complex order parameter fields^ the actual crit- 
ical behavior at low N and in the strong coupling limit 
may be compatible with a second-order phase transition. 
For a GL model with a single complex order field, this 
expectation is confirmed both by a duality and MC anal- 
ysis of the modeh 22 ' 23 ' 24 ' 25 ' 26 ' 27 ' 28 There are RG stud- 
ies where low- TV critical points were found ; 29 ' 30 ! 31 ' 32 but 
they all have some kind of problem, being either not well 
controlled or involving ad hoc assumptions. The two- 



loop perturbation in terms of the e-expansion should be 
in principle a well controlled approach. 33 The resumma- 
tion of the two-loop result^ , led to a critical point for 
N = 1. However, the unresummed result, though cor- 
rect, contains unacceptable pathologies, like for example 
a critical exponent v for 182.9 < N < 200 larger than the 
N — > 00 result at fixed dimension, and the absence of a 
fixed point for the gauge coupling if N < 18e. In view of 
these pathologies, the resummed two-loop result cannot 
be completely trusted. The ideal approach would be to 
obtain a resummed three- loop result. Unfortunately, to 
this order the only available RG function is the (3 func- 
tion for the gauge coupling. 34 Anyway, the important ar- 
gument from Ref. [Tt] to be retained here is the following. 
The existence of a critical value of N in the weak cou- 
pling analysis in d = 4 — e spacetime dimensions, above 
which a critical point exists, is thought as an indication 
that in the strong-coupling limit a quantum critical point 
at lower values of N may exist. This is of course a conjec- 
ture that must be tested further. Note, however, that the 
easy-plane case, also studied in Ref. Qj], when generalized 
to a theory with a global O(N) x O(N) symmetry, does 
not have any fixed points at nonzero gauge coupling for 
all values of N, i.e., this theory does not have a critical 
value of N. 

The behavior of the e-expansion for the GL model with 
a single complex order field improves considerably if the 
theory is coupled via the gauge field to Nf Dirac fermions 
species^ Interestingly, it is not necessary to have a large 
value of Nf to obtain an infrared stable fixed point. In 
the context of the present paper, it should be noted that 
the SU(2) model for deconfined spinons is exactly the 
same as a GL model with two complex order fields. 

In this paper we will consider the Lagrangian: 



C = Cb + Cf, 



(1) 



where 



1 N b 
£b = ■^-2(e fl uxd v Ax) 2 + ^2 \(dn - iA^)z a \ 



and 



N b / N b \ 



N f 



(2) 



(3) 



The Lagrangian Cb is precisely the model for deconfined 
spinons proposed in Ref. 0. In this context, we are 
primarily interested in the case with Nb = 2. The La- 
grangian Cb is supposed to govern the universality class 
of the phase transition between a Neel state and a VBS. 
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The Lagrangian Cf contains Nt species of four- 
component Dirac fermions in 2 + 1 Euclidean 
dimensions.— The representation of 7 matrices in 
this case is given by 
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(4) 



where o\, a 2 , and 03 are the Pauli matrices. 

The Lagrangian Cf is typical of a so called algebraic 
quantum liquid. There are three such quantum liquids, 
the algebraic spin liquid (ASL),— the algebraic Fermi 
liquid (AFL)j2£ and the recently introduced algebraic 
charge liquid (ACL)4^ Let us briefly mention the dif- 
ferences between these algebraic quantum liquids. 

The ASL emerges out of the so called staggered flux 
phased in the large Nf limit of a SU (Nf) Heisenberg an- 
tiferromagnet. In this case the spin operators are writ- 
ten in terms of fermion bilinears such that the Heisen- 
berg model reads H = -(J '/N f )Y, {itj) fiJjafjpfip, 

with the constraint f} a fia — Nf/2. A lattice gauge 
field arises from the phase fluctuations of the Hubbard- 
Stratonovich link field Xij — (fl a fja)- This gives rise to 
a compact U(l) gauge theory. Linearizing around the 
nodes ±(tt/2, 7r/2) of the quasi-paricle spectrum — 
2| Xo I \/ cos 2 k x + cos 2 k y and taking the continuum limit 
leads to the Lagrangian ([3|). For large enough Nf the 
fermionic spinons deconfin e 42 i 43 ' 44 and we obtain a Mott 
insulating state having no broken symmetries, i.e., a spin 
liquid. For low enough Nf the spin liquid is no longer 
stable 44 and the chiral symmetry is probably broken due 
to spinon confinement. However, in this theory chiral 
symmetry breaking (CSB) takes place even in the non- 
compact casei 36 ' 45 This CSB is associated to the devel- 
opment of an insulating antiferromagnetic stated 

The AFL, though having essentially the same La- 
grangian Cf, is physically very different from the ASL. 
Here the gauge field is not compact and has a completely 
different origin. 38 i 39 While the ASL is associated to a 
Mott insulator, the AFL originates from a <i-wave super- 
conductor. The Dirac fermions are obtained by identi- 
fying the low-energy quasi-particles on the nodes of the 
d-wave gap. The gauge field follows from the coupling 
of these quasi-particles to vortices. The way this is done 
is very subtle and involves the manipulation of singular 
gauge transformations. 3 - The fermions in the AFL, just 
like the ones of the ASL, carry no charge. Interestingly, 
the number Nf in the AFL is related to the number of 
layers of the system. For example, a single-layer system 
has Nf — 2, while a bilayer system has Nf = 4. Note 
that in this case there is no need for other gauge fields in 
the corresponding layers, since the gauge field arises from 
a vortex-antivortex excitation. Thus, at large distances 



a vortex in one layer is always connected with a vortex in 
the second layer, the same happening for the antivortices. 
The end result is that Dirac fermions in different layers 
couple to the same gauge field4£ Of course, similarly to 
the ASL, CSB also occurs here for small enough A/ t 39 i 48 
and is once more associated to antiferromagnetism. 

In the ACL a fractional fermionic particle with charge 
e and no spin, a holon, couples to bosonic spinons via 
a gauge field4°- Therefore, it is necessarily related to the 
concept of DQC and in this case the whole Lagrangian (fl]) 
has to be considered. Note that in contrast with the ASL 
and AFL, the Dirac fermions in the ACL carry charge. 
This new state of matter bears some resemblances with 
earlier ideas on spin-charge separation in the cupratesjS 
An important difference is that in earlier theories super- 
conductivity is obtained by doping a spin liquid, while 
in the ACL superconductivity arises by doping an anti- 
ferromagnetic state near the deconfined quantum critical 
point. 49 Thus, we can either dope a Neel or a VBS state. 
Like in the two other algebraic quantum liquids, here 
the chiral symmetry can also be broken. However, since 
the Dirac fermions now carry charge, we have a situation 
where charge density waves develop as a result of CSB. 
It is important to note that the coupling to bosons leads 
to a reduction of the value of Nf below which the chiral 
symmetry breaks 4£ 

From the above discussion we see that among the three 
algebraic quantum liquids, only the ACL cannot exist 
without the coupling to bosons. The results of this pa- 
per will concern mainly this case. Quantum criticality 
in U(l) gauge theories with both bosonic and fermionic 
matter has been also considered in a recent paper^ 
There the authors considered a CP^ 6-1 model coupled 
to Nf Dirac fermion species and analyzed the model for 
large Nf, and Nf, while keeping Nf/Nb arbitrary. Here 
we have followed Ref. H and softened the CP^ 6-1 con- 
straint. This leads to a theory that is tractable in a 
RG framework in d = 4 — e spacetime dimensions. The 
main advantage of this approach is that we will be able 
to work with a fixed Nf, = 2 and Nf does not need to 
be large to obtain a quantum critical point. Part of our 
results follow directly from Ref. [HI, the main difference 
being that here we are interested in Ah = 2 rather than 
Ah = l. 52 However, in this paper we will improve sub- 
stantially upon the previous analysis and present many 
new results which are relevant for the quantum critical 
behavior of the ACL. One of the main new results of this 
paper will be the calculation of i]n for Nb — 2 and Nf > 4 
in terms of the crossover exponent tp. As shown in Ref. 
UtI the critical exponent 77^ is related to the crossover 
exponent by the formula 



r] N = d + 2(l - ip/v), 



(5) 



where v is the correlation length exponent. The crossover 
exponent is related to a mass anisotropy of the La- 
grangian. Although the Lagrangian Cb does not have 
any mass anisotropy, the correlation function (z*(x) • 
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z(0) z(x) • z*(0)) is calculated through the insertion of 
the operator z* a zp. For a ^ f3 this insertion allows us to 
calculate (/?.— 

The plan of the paper is as follows. In order to ex- 
plain why a second-order phase transition is easier to 
obtain when Dirac fermions are included, in Section II 
we will calculate the effective potential associated to the 
Lagrangian (fTJ). There we compare the cases with and 
without fermions. In the absence of fermions the one- 
loop effective potential typically describes a first-order 
phase transition. When fermions are included the situa- 
tion changes and we can show that for small Nf a first- 
order transition takes place, while for larger values of Nf 
the effective potential features a second-order phase tran- 
sition. In Section III we proceed with the RG analysis 
and the calculation of the critical exponents. In Section 
IV we briefly discuss the effect of CSB and show that it 
does not affect the quantum critical regime analyzed in 
Section III. We also compute the anomalous dimension of 
the chiral holon susceptibility. Section V concludes the 
paper. An Appendix sketches the details of the calcula- 
tions of the anomalous dimension tjn- 



where A is a ultraviolet cutoff and we have neglected 
terms which are independent from (j>. The above effec- 
tive potential can be written in terms of renormalized 
quantities m 2 and u defined by the normalization con- 
ditions J7" ff (0) = m 2 and 1/^(0) —_u, where the primes 
denote derivatives with respect to <fi. This leads to 



and 



rn 



u = u 
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2^"' 



3^ f 
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(9) 



(10) 



where in Eq. (|10p we have replaced in the second term 
m 2 , by m 2 , since the error involved in this replacement 
contributes to an order higher than the one being calcu- 
lated here. Thus, the effective potential becomes 



m2 A? , U ° 14 1 ( 2 . u 

4) H 7<t> m H c 

2 V 4! v 12tt V 2 s 
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(11) 



II. EFFECTIVE POTENTIALS AND THE 
ORDER OF THE PHASE TRANSITION 

A. A simple example 



The above potential is typical of a system exhibiting a 
second-order phase transition. The order parameter is 
obtained by extremizing the effective potential. It is 
given by 



In order to illustrate the subtlety involving the deter- 
mination of the order of the phase transition in models for 
deconfined spinons, we will consider the following simple 
Lagrangian for an interacting scalar theory in d = 2 + 1 
Euclidean dimensions: 



£ = 2^ 



,2 , U 



(6) 



At the mean-field level the above model obviously ex- 
hibits a second-order phase transition with a critical 
point at tuq — 0. Let us calculate the effective poten- 
tial for the above model at one-loop order. This is more 
easily done by writing </> = </> + 5(p, where <j> is a constant 
background field, and integrating out the quadratic fluc- 
tuations in 8<j) while disregarding the higher order ones. 
This calculation is very simple, since it just involves a 
Gaussian integration in 6<fi. We obtain, 



= ± 



V3_ 

8ituq 



3uZ 



1287r 2 m 2 w 



9ug 



1/2 



5127r 2 m 2 UQ 
(12) 

and we see that the critical point is given by m 2 = 0. 
Looking at Eq. (TTU)) we may think that u will get large 
and negative as the critical point is approached, such that 
the effective potential would become unstable against a 
<fi 6 interaction. However, this is not the case, since up to 
one-loop accuracy we can write the following RG equa- 
tion for the dimensionless coupling u — u/(16irm), 



du „ 2 
m— — = — u + 3u . 
dm 



(13) 



We see that as m — > the above (3 function vanishes at 
the infrared stable fixed point = 1/3. Thus, u actually 
does not diverge as the critical point is approached. 



+ ^ [Trln(-d 2 + m 2 + 2f0 2 ) -Trln(-9 2 )] , (7) 

where V is the (infinite) volume. Explicit evaluation of 
the tracelog term yields 

TT fl\ 1 ( 2 . Au o\ 72/0 74 1 ( 2 . U 7 2 \ 3 

(8) 



B. Effective potential for the SU{2) 
ant ifer romagnet 

Let us consider now the one-loop effective potential 
for the Lagrangian ([2]). Like in the example of a single 
scalar field theory, all we have to do is to integrate out 
the Gaussian fluctuations. Note that the Lagrangian is 
already quadratic in the gauge field, so that the latter can 
' be immediately integrated out exactly. We will consider 
a nonzero background for the fields z\ and Z2, leaving 
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a vanishing background for the remaining iV& — 2 fields. background spin orientation field n = z a cr a pzp. The 
Thus, we will have Z\ = z \ + <5zi, z 2 = z 2 + < 5- z 2 J an d result is 
z Q = 5z a for a > 3. Thus, we have for Nf, = 2 the 



U^( Zl ,Z2) = m 2 (N 2 + N 2 ) + l -f(\^ + |z 2 | 2 ) 2 - # (N 2 + |z 2 | 2 ) 3/2 
^ { [m 2 + MQ (3|zi| 2 + N 2 )] 3/2 + K + Uo (N 2 + 3|z 2 | 2 )] 3/2 + 2(7V b - 1) [m 2 + Uo (N 2 + N 2 )] 3/2 } 



where m 2 is given by 



m = r + 



-no 



4e 2 



A 



(14) 



(15) 



The third term in Eq. (|14p arises from the integration over the gauge field and is reminiscent from the fluctuation- 
corrected mean-field theory in the GL superconductor j 21 1 54 Due to this term the symmetry will be broken for a 
in 2 > 0, leading in this way to a weak first-order phase transition. 2 — Note that for eo = we have a second-order 
phase transition in the universality class of a 0{2N\ } ) symmetric classical magnetic system in three dimensions. 

C. Effective potential for the algebraic charge liquid 

Finally, let us consider the the one-loop effective potential for the the Lagrangian (fTJ). First we note that the 
one-loop contribution following from integrating out the Dirac fermions is given by 



(16) 



where = d^A v — dyA^ (note that the above expression is multiplied by a factor two in the case of two-component 
Dirac fermions). The above expression follows simply by calculating the one- loop vacuum polarization in 2 + 1 
dimensions. The background fields are chosen in the same way as in the previous Subsection. After integrating out 
the quadratic gauge fluctuations for a fixed spinon background, we obtain the correction 



1 



6Cf = -^Trln 



-9 2 + ^Pv / ^+2e 2 (|z 1 | 2 + |z 2 | 2 ) 



Trln(-<9 2 ) 

After evaluating the above tracelog and integrating out the spinon Gaussian fluctuations, we obtain 



(17) 



U^^) = m 2 (N 2 + N 2 ) + f (| 2 "i| 2 + N 2 ) 2 + I :,.| 



(iV /e 2 /8) 3 - (5iY /eo 74)(N 2 + N 2 ) + (128e 2 /7V / )(|z 1 | 2 + \z 2 ^ 



384 



In 



2e 2 (N 2 + |z 2 | 2 ) 



A 2 



12tp 
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"f e 
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In 



l-Jl 



#f,(N 2 + N 2 ) 
f 



512 
"f e 



lk { [m ' + u ° (3|z_l12 + |2 " 2|2)]3/2 + [m2 + u ° (|fl ' 2 + 3 I Z " 2 ' 2 )] 3/2 + W - 1) K + «o (N 2 + N 2 )] 3/2 } , 

(18) 

where m 2 is still given by Eq. (fT5|) . Note that the annoying third term of Eq. (|14[) disappeared once Dirac fermions 
were included. For Nf — ► the above equation reduces to Eq. (jT4j) , as it should. Thus, we can expect that for 
Nf sufficiently small a first-order phase transition occurs. On the other hand, for Nf above some critical value a 
second-order phase transition should take place. The leading small Nf correction shifts the mass term and add a 
logarithmic correction, resulting in 
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(19) 



The phase transition described by the above effective potential is certainly a first-order one. 

To see how the second-order phase transition takes place, we perform a large Nf expansion with the limit Nf — > oo 
being taken with Nfe^ kept fixed. Up to order l/Nf we have 



{ [m 2 + Uo (3|zi| 2 + N 2 )] 3/2 + [m 2 + uo ( | ^ | 2 + 3|z 2 | 2 )] 3/2 + 2(N b 1) [m 2 + (N 2 + %?)f 2 } 



(20) 



where 
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The (|zi| 2 + |z2| 2 ) 2 -term is reminiscent of the Landau 
expansion of the dual formulation of the GL model 
of a superconductor ) 23 ' 51 where a tricritical point was 
shown to exist. By defining a "Ginzburg parameter" 
k 2 = uo/(2e 2 l ), we see that the effective potential (|2"0")) 
is stable only for k 2 > 8/(Nfir 2 ), in which case a second- 
order phase transition occurs. 

The precise characterization of the change of behavior 
as a function of Nf will be studied in the next Section 
by means of the RG. 



III. QUANTUM CRITICAL BEHAVIOR OF 
THE SU(2) ALGEBRAIC CHARGE LIQUID 

The P functions at one-loop order in d = 4 — e dimen- 
sions are easily calculated using standard methodsj 11 ' 55 
The dimensionless couplings are given as a function of 
the momentum scale fj, by / = ZA^~ e &^ and g — 
Z 2 /j l ~' i uo/Z Ul where Za and Z z are the wave function 
renormalizations for the gauge field and spinon field, re- 
spectively. Z u is obtained from the vertex function as- 
sociated to spinon self-interaction. The renormalization 
factor Za is obtained directly from the one-loop vac- 
uum polarization. The f3 functions /?/ = [idf/d/i and 
j3 g = jidg/d^i are given by 



Pf = -ef + 



4N f +N b f2 ^ 



and 



1% 



-eg-6fg+(N b + A)g 2 + 6f, 



(22) 



(23) 



where / and g were rescaled as / — > 87r 2 / and g — > 8Tr 2 g 
in order to remove unnecessary geometrical factors. Note 
that at one-loop the number of fermion components Nf 



(21) appears only in Eq. 



Of course, for N 



J 



the 



above result coincides with the one for the Ginzburg- 
Landau superconductor i 21 i 30 i 33 

For Nb = 2, and provided Nf > 4, there are two fixed 
points (g±, /*) at nonzero gauge coupling whose coordi- 
nates are given by 



/* 



3e 



2(2JV) + 1)' 



(24) 



and 



5 + N f ± 



9± = 



N f 



10Nf - 56 



6(2JV/ + 1) 



(25) 



From the above fixed points, only (g+, /*) is infrared sta- 
ble, thus corresponding to the quantum critical point 
where the critical exponents will be calculated. A 
schematic flow diagram is shown in Fig. 1. We see that 
the effect of Dirac fermions on the phase structure is quite 
remarkable, since neither N b nor Nf need to be large in 
order to attain criticality. This behavior is in strong con- 
trast with the one obtained in the limit Nf — > 0, where 
quantum criticality occurs only for N b > 182. 9f^ 

The case Nt = 4 exhibits a peculiar critical behavior, 
as the calculations of the critical exponents will show. In 
fact, for Nf = 4 the fixed points (g-,f*) and (g+,f*) 
coincide. 

The critical exponent v is obtained through the cal- 
culation of the insertion of the operator ^ Q |z a | 2 in the 
correlation function (z a (x)z^(y))i 11 ^ 55 This leads to new 
singularities, so that we have to introduce a new renor- 
malization constant, Z 2 , associated with this operator 
insertion. In pure scalar theories Z 2 is very simply cal- 
culated by differentiating the two-point function at zero 
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0(4) 
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FIG. 1: Schematic flow diagram associated with the j3 func- 
tions and (J23J) for jV 6 = 2 and Nf > 4. In the figure 
P± = (g± , /* ) . Note that only P+ is infrared stable, thus 
corresponding to a quantum critical point. For a vanishing 
gauge coupling we have a second-order phase transition gov- 
erned by a 0(4) Wilson-Fisher fixed point. This fixed point 
is unstable for / =fc 0. 



Q(x) = 2{z*(x) ■ z(0) z(x) ■ z*(0)) - — (|z(*)| 2 |z(0)| 2 ). 

(30) 

For Nb — 2 we have simply 



Q{x) = 2(z* 1 (x)z 1 (0)z 2 (x)z* 2 (0)) +2(z 1 (x)z* 1 (0)z* 2 (x)z 2 {0)) 
~ 2<|z 1 ( a ;)| 2 |z 2 (0)| 2 ) + <|z 1 (x)| 2 |z 1 (0)| 2 ) 
+ (|z 2 (x)| 2 |z 2 (0)| 2 ). (31) 

The critical exponent 77 at is defined by the behavior of 
Q{x) at the quantum critical point: 



Q{x) 



1 



Md— 2+rjjv 



(32) 



where t]n is given in terms of the crossover exponent if 
in Eq. 

We also note the following scaling behavior at the 
quantum critical point: 



momentum with respect to the mass. However, in gen- 
eral is better to perform the renormalization at nonzero 
momenta in order to avoid certain infrared divergences, 
especially when the scalar fields are coupled to a gauge 
field, like the case studied here. The standard way to 
calculate v is to compute the following RG function at 
the infrared stable fixed point: 



d\n{Z 2 /Z z ) 

72 = 75 : 

an 



(26) 



which leads to 



»72, 



(27) 



where 772 is the value of 72 at the infrared stable fixed 
point. At one-loop order, we have 



72 = 3/ - (N b + l)g. 



(28) 



By inserting /* and g + for Nb = 2 in the above equations 
and expanding up to first order in e, we obtain 



Nt 



Nj + 10N f - 56 



4(l + 2JVy) 



0(e 2 ). (29) 



The critical exponent 77 at is also derived through a 
quadratic operator insertion. However, in this case there 
is a mass anisotropy involved. For a general Nj, the cor- 
relation function Q{x) = (n(x) • n(0)) is given in terms of 
spinon fields as 



X)M*)| 2 X>(o)| s 



la.ld-2+774 ' 



(33) 



where 



774 =d + 2(l- \jv). 



(34) 



Thus the scaling dimension of the operator ^ Q |z q (;e)| 2 
is given by 



dim 



l/v. 



(35) 



That 774 depends only on v is very easy to see, since it 
follows from the correlation function between operators 
So I ^q, 1 2 , whose singular behavior leads to the renormal- 
ization constant Z 2 and consequently to v. From the 
expression (I29p for the critical exponent z/, we obtain 



774 = 2 



N) 



107V/ - 56 



1 + 27V/ 



e + 0(e 2 ) (36) 



The calculation of 77^ depends on the calculation of the 
crossover exponent if. The calculation is facilitated by 
observing that the scaling dimension of each component 
of n is given by d — fjv. In particular, we have 



dirndl 2 - |z 2 | 2 ] = d-ip/v. 



(37) 



Thus, ip can be calculated from the insertions of |zi| 2 and 
\z 2 \ 2 in the correlation function Gu(x) = (zi(x)zl(0))£& 
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TABLE I: Values of the critical exponents for Nt — 2 and 
different values of Nf > 4. We have set e = 1. 









77 i 

f/4 




4 


1 19 

1/Z 


1/ .3 


I 


1^/19 


K 



fl 7Q 




u.ooyzy 


1 071 1 7 
l.Z / 11 / 


6 


0.66009 


0.75191 


0.33468 


1.3245 


7 


0.68229 


0.84305 


0.3417 


1.35381 


8 


0.69678 


0.90943 


0.36696 


1.37208 


9 


0.70689 


0.96007 


0.39749 


1.38429 


10 


5/7 


1 


3/7 


39/28 


11 


0.71988 


1.0323 


0.45837 


1.39907 


20 


0.73978 


1.17336 


0.64274 


1.41792 


50 


0.74817 


1.27148 


0.83646 


1.42103 



This calculation is sketched in Appendix A. By denoting 
those insertions by Gu a , where a = 1,2, we can ob- 
tain the crossover exponent from the renormalization of 
Gii,l{x) — Gn t 2{x). This procedure leads to a new renor- 
malization constant Z' 2 and, analogously to Eq. (|26[) . we 
can define the RG function 



7 2 = 



dln(Z' 2 /Z z 



dp 



which at the infrared stable fixed point yields 



v 



(38) 



(39) 



By comparing the above equation with Eq. (|5|) we see 
that t]n — d — 2 — 2772 • At one-loop order we have 



72 = 3/ - 9, 



and therefore 



(40) 



N 2 + 10N f - 56 - 5AT/ - 16 
2 + — rrr- e+0(e 2 ), (41) 



3(1 + 2N f ) 
while the crossover exponent is given by 



(N f + ^JN 2 } + WNf - 56) 
12(1 + 2^) 



56-11 



0(e 2 ). 



(42) 

The critical exponents for three spacetime dimensions 
(e = 1) are shown for several values of Nf at Table I. As 
already mentioned, the case Nf = 4 is peculiar. Indeed, 
for Nf = 4 the critical exponents are rational and, more- 
over, it corresponds to the only case where r/^ < rj4,. Note 
that for Nf — 10 the exponents are also rational, but in 
this case 77 at > ?74, just like for all the other values of Nf. 
Interestingly, we have that ^jv > 1 for Nf > 10. This 
result is consistent with the analysis at large Nf and 
with Nb/Nf arbitrary. 50 Indeed, in that case we find that 
77.7V > 1 if Nf /Nb > 3 (if four-component Dirac fermions 
are used, the latter inequality becomes Nf/Nj, > 3/2). 

IV. CHIRAL SYMMETRY BREAKING AND 
CHIRAL SUSCEPTIBILITY 

We have already mentioned in the Introduction that 
CSB constitutes an important physical aspect of both 
algebraic Fermi and spin liquids. In this Section we want 
to investigate the effect of CSB on the ACL^ This is 
typically a phenomenon that takes place at low energies. 
In the case of AFL and ASL it occurs for \p\ < N f el^^ 
For the ACL this corresponds roughly to \p\ <C a = (Nf + 
Nb/2)eQ. In this regime the gauge field propagator is 
dominated by the one-loop vacuum polarization and is 
given by 



DM 



16 



(2N f + N b )e 2 \p\ 



p 2 



(43) 



where we have used the Landau gauge and assumed a 
three-dimensional Euclidean spacetime. Following the 
standard literature ; 36 ^ 5 we consider the one-loop fermion 
Green function with a dressed fermionic propagator and 
where the above gauge field propagator has been in- 
serted. The inverse propagator is given by G^ ! 1 (p) — 
ip'Z(p) + 'E(p) and the corresponding self-consistent equa- 
tion reads 



r -i,s = 16 f f d 3 k 1 ^(k~p) + i^-j>)Z{k-p)] 1 ^ _ f d 3 k t[E(k-p) + itf-t)Z(k-p)]fk \ 

W> W ^A^ + A^y (27r) 3 [Z 2 (Jfe-p)(/c-p) 2 + E 2 (fc-p)]|fc| J (2tt)3 \Z 2 (k - p)(k - p) 2 + Z 2 (k - p)]\k\ s j 

(44) 

Note that the only difference with respect to the standard analysis is the coefficient in front of the curly brackets, 
which gets modified due to the spinon loop contribution to the vacuum polarization. We will set Nb = 2, which is the 
case of physical interest to us. Futhermore, in order to simplify the analysis we will assume that Z(p) ss h 45 i 57 After 
standard manipulations involving the algebra of the gamma matrices, we obtain the self-consistent equation for the 
self-energy: 



9 



S(p) 



32 



2N f 



Integrating over the angles, we obtain 



d 3 k 



E(fe) 



(2tt) 3 [k 2 + E 2 {k)]\k+p\ 



(45) 



(2N f + l)ir 2 \p\ J 



16 

(2JV/ + 1)tt 2 |p| 



d\k 



|*|E(*)(|A| + |p|-|*-p|) 



fc 2 + E 2 (fc) 



k 2 + E 2 (fc) 
I 



fc|E(fc) 



bl 



E 2 (fc) 



(46) 



which can easily be converted to the differential equation 



d 
dp 



dp 



16 



p 2 S(p) 



7r 2 (2iV/ + l)p 2 + E 2 (p)' 



(47) 



The above equation can be linearized in the regime 
\p\ ^S> E(p) where it can be solved. We will omit the de- 
tails here, since the analysis is well known.— The solution 
here is only slightly changed due to the spinon degrees 
of freedom. The result is a dynamically generated mass 
gap E(0) ~ (ipip) of the form 



E(0) = a exp 



2tt 



64 

7T 2 (2JV/+1) 



1 



(48) 



This mass generation is associated to the development 
of a chiral condensate, which signals the appearence of 
holon charge density waves in the system. The above 
mass gap vanishes for Nf > N c h, where 



which defines the critical exponent 774 . If fluctuations are 
ignored, we can write simply 



G x (x) 



1 



\ x \2(,d-l)' 



(52) 



which leads to 774 = d. 

We can obtain the anomalous dimension 774 in 2 + 1 
dimensions for large Nf and Nb, with Nf/Nf, arbitrary, 
directly from the result of Ref. |60| by the simple replace- 
ment 8/N -> 8/(N f + N b /2). The result agrees of course 
with the one obtained in Ref. IBol . up to a trivial differ- 
ence related to the fact that the latter reference works 
with two-component Dirac fermions. 

Let us compute now 774 for the SU(2) ACL in first- 
order in e. Similarly to the case of four-spinon correlation 
functions, the calculation can be done by computing the 
insertion of the operator ipip in the correlation function 
G^(x). This amounts in computing the scaling dimension 
of the operator ipip, i.e., 



N, 



ch 



64 
^2 



1 



2.74, 



(49) 



which is smaller than the corresponding value in ab- 
sence of spinonSj2§ Thus, CSB seems to occur only in the 
regime where no quantum critical point is found. There- 
fore, by accepting the calculations and underlying ap- 
proximations employed so far in this paper, it seems that 
CSB does not spoil the quantum critical behavior. 

Another correlation function of interest is the chiral 
susceptibility 



G x {x) = (iP(x)iP(x)iP(0)tP(0)), 



(50) 



whose scaling behavior for a large number of fermion 
components and 2 + 1 dimensions has been calculated 
in Refs. [5^ and [60l . At the quantum critical point the 
chiral susceptibility behaves like 



G x (x) 



1 



Uld-2+rji ' 



(51) 



dim[^^] 



1 



(53) 



where 772 is the fixed point value of the RG function re- 
lated to mass renormalization. Thus, the easiest way to 
compute the above scaling dimension is by adding a bare 
mass term in the Lagrangian Cf. The calculation of the 
needed renormalization constant is standard and can be 
found in quantum field theory textbooks. 55 We obtain at 
one-loop the result t) 2 = 3/*. From Eqs. |5T|) and (|53|) . 
we obtain 



and therefore, 



774 = 2 



fji = d- 2fj 2 , 



2-^±e + 0(e>) 



2AT/ + 1 



(54) 



(55) 



In particular, for Nf = 4 and e = 1 we obtain 774 = 2. In 
this case we have that the chiral susceptibility behaves 
at the quantum critical point as 
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(56) 



which leads to a logarithmic behavior in momentum 
space. 

V. CONCLUSION AND DISCUSSION 

In this paper we have studied the quantum critical be- 
havior of a field theory associated to the algebraic charge 
liquid (abreviated as ACL throughout this paper), which 
is a new type of quantum liquid proposed recently in 
Ref. HO- Our study was made in the framework of the 
e-expansion. To lowest order it leads to a quantum crit- 
ical point provided the number of fermion components 
Nt > 4. Since we are considering a representation using 
four-component spinors, the physically relevant number 
of fermions is Nf = 2M- Thus, in the framework of the 
present study, we are unable to find a deconfined quan- 
tum critical point for Nf = 2. However, it is likely that 
an improved approximation will show that the theory 
is also critical for Nf = 2. The point is that the e- 
expansion for the Lagrangian should, in a nonpertur- 
bative framework, exhibit an infrared stable fixed point 
for all values of JVj, > 1. The e-expansion cannot cap- 
ture such a strong-coupling regime. At the moment there 
is no reliable nonperturbative scheme for this problem. 
As discussed in the introduction, a resummed three-loop 
e-expansion would provide the desirable controlled ap- 
proximation. The inclusion of Dirac fermions make the 
e-expansion more reliable, since they essentially make the 
gauge coupling weaker. It is remarkable that in such a 
case it is not any longer necessary to have a large number 
of spinons to obtain a quantum critical point. Interest- 
ingly, the number of Dirac fermions does not need to be 
large cither. Thus, we were able to explicitly exhibit a 
deconfined quantum critical point for an SU (2) theory of 
bosonic spinons. 

In order to convincingly show how the Dirac fermions 
make the full theory more well behaved, let us have a 
look at the two-loop j3 function for the gauge coupling. 
This result can be easily derived by using the earlier two- 
loop result for the Lagrangian combined with the 
higher order result for the Lagrangian ([3]). 62,63 The result 
is 



Pf = ~*f + Nb \ ANf f 2 + 2(N b + 2N f )f\ (57) 



that we still have an infrared stable fixed point only for 
N b > 182.9. Thus, besides the two-loop result for N f = 
being unable to produce a quantum critical point for low 
values of 2V&, it introduces a new difficulty, since a non- 
trivial fixed point for / does not any longer exists for all 
Nb > 1. However, the situation improves considerably 
if Nf 0. Indeed, in this case we have that /* > 
for N b > 18e - AN f . Thus, for N b = 2 we should have 
Nf > (9e — l)/2, which after setting e = 1 becomes 
Nf > 4. While this clearly shows that the introduction 
of Dirac fermions improves the behavior of the SU (2) an- 
tiferromagnet, it is a little bit frustrating to see that at 
two-loop no quantum critical point for Nf = 4 is found 
anymore. However, we should keep in mind that the the- 
ory in absence of fermions is not being considered in a 
strong-coupling regime. Once this is done in some way, 
a critical regime even for values as low as Nf = 2 should 
be expected to occur. Note that even if such a result is 
achieved, we may have to face another problem, namely, 
CSB. The calculation in Sect. IV indicates that CSB 
leads to gapped fermions for Nf = 2. If this number is 
really correct, then there will be no deconfined quantum 
critical point anyway. Physically this means that CDW 
will screen the interspinon interaction in such a way as to 
make the spinon correlation length finite, and the only 
possible phase transition, if any, would be a first-order 
one. Adopting a more optimistic point of view, let us 
mention that there is strong evidence that the values of 
Ach calculated using methods like the one we have used 
in Sect. IV tend to be overestimatedj 64 i 65 ' 66 Since the 
spinons contribute further to reduce the value of A c h, we 
may hope that the quantum critical point for Nf = 2 will 
survive. 

Larger values of Nf are physically relevant for multi- 
layer systems. As discussed in the Introduction, this is 
precisely the case with the AFL.— However, we must be 
aware of the fact that in the case of the ACL additional 
gauge fields may play a role, so that the analysis made 
here would have to be modified accordingly. 

The main achievement of this paper was the computa- 
tion of the anomalous dimensions of composite operators 
associated to the quantum critical behavior of the La- 
grangian {1} . Particularly important was the calculation 
of the anomalous dimension rj^ of the Neel field. The 
results presented here constitute an explicit and well con- 
trolled example of deconfined quantum critical point for 
a system having SU (2) invariance. The computed values 
of t]n are substantially larger than the ones that would 
follow from a LGW analysis, reflecting a typical feature 
of the DQC scenario. 



which has a nontrivial fixed point at 



/* = 



3e 



N h + AN S 



1 - 



18 



N b + ANf 



e + 0{e 2 ) 



(58) 



For Nf — we obtain that /* > only for N b > 18e. Fur- 
thermore, in this limit the two-loop result for (3 g shows 
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FIG. 2: One-loop Feynman diagrams representing the inser- 
tion of the operator \zi\ 2 (dashed line) in the correlation func- 
tion (zi(x)zi(0)). Here the lines labeled 1 or 2 are associated 
to the fields Z\ and Z2 respectively. The wavy line represents 
the gauge field propagator. 




1 



1 



FIG. 3: For the insertion of the operator |z2| 2 (dashed line) in 
the correlation function (zi(x)z*(0)) only one diagram con- 
tributes. 



APPENDIX A: CALCULATION OF m 

In this appendix we sketch the calculation of r/N based 
on the insertion of the operators \zi\ 2 and \z2\ 2 in the cor- 
relation function Gn(x) = (zi (x)z* (0)). The calculation 
follows the method of Ref. [53J, i.e., it is based on dimen- 
sional regularization in the minimal subtraction scheme. 

In Fig. 2 the one-loop Feynman diagrams associated 
to the insertion of \zi\ 2 are shown. For the insertion of 
| ^2 1 2 , on the other hand, only one diagram is produced 
at one-loop (Fig. 3). 

The renormalization constant Z' 2 is determined by the 
normalization condition in momentum space 



z 2 [Gn,i(pi,Pr,Pa.) - G 11 1 2 (p 1 ,p2;p 3 )] \ S p = 1, (Al) 

where p$ = —(P1+P2) and SP denotes the symmetry 
point defined by 



Pi ■ Pj 



(A2) 



We will also need to compute Z z , which can be deter- 
mined from the normalization condition 



Z 



dp 2 



= 1. 



(A3) 



A straightforward way to understand the insertion of 



the operators 
of the form 



zi\ and I ^2 1 is by adding a source term 



^source = ri(x)|zi(a;)| 2 + r 2 {x)\z 2 (x)Y 
to the Lagrangian. Thus, 



(A4) 



in perturbation theory. Thus, the ratio Z 2 /Z z is gauge 
independent, which should be expected physically, since 
it leads to the anomalous dimension of the gauge invari- 
ant operator |zi| 2 — | ^2 1 2 - A similar argument holds for 
Z 2 /Z z , which leads to the critical exponent v. Since 
Z' 2 jZ z is gauge independent, any gauge can be fixed to 
calculate Z z and Z' 2 . Here we are going to use the Feyn- 
man gauge, so that the gauge field propagator is given in 
imaginary time by 



(A6) 



At the quantum critical point and at one- loop order, 
we have 



G^(p) = P 2 



2 2 f d d k fipp - kp)(2p v - k v ) 



= p* 1 - 



' (2ir) d 
el\p\ 



4ir 2 c 



(p — k) 2 
finite, 



(A7) 



where we have made use of the integral 



d d k 



T(2-d/2)T 2 (d/2-l)^ d _ 4 



(2w) d (p-k) 2 k 2 



(47r) d / 2 r(d - 2) 



(A8) 

and the usual properties of dimensional 
regularization! 11 ! 55 Next we use Eq. (|A3|) to obtain 



Z z = 1 



47r 2 e ' 



(A9) 



which can be rewritten in terms of the renormalized di- 
mensionless coupling up to one-loop accuracy as 



Z z = 1 



2/ 



(A10) 



Gn. a {x\x') = 



5r a (x') 



(A5) 



Note that both Z z and Z' 2 are gauge dependent. How- 
ever, this gauge dependence can be calculated exactly 
from the Ward identities^ 5 - It can be shown that Z z and 
Z' 2 have both the same gauge dependence to all orders 



The calculation of Z 2 is made similarly. In this case 
standard manipulations in Eq. (|Alj) leads once more to 
integrals of the form (|A8[) . An integral of the form 



d d k 



(27r) d ( Pl -k) 2 (p 2 + k) 2 k 2 



(AH) 
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also appears in Eq. (|A1[) . However, there is no need 
to worry with the above integral, since it does not have 
poles at e = 0. At the end we obtain the result 



The renormalization constant Z' 2 follows from 



Therefore, 



z' 2 = i + g --L 

2 e e 



Z> 



9 3/ 



(A12) 



(A13) 



[Gu4(0,p; -P) - G^] 2 (0,p; -p)] \ ri=r2= ^ = Z' 2 \ 

(A17) 

and the result (|A12[) is obtained once more. 

In order to obtain 77 jy, we first note that the renor- 
malized direction field, which is a composite operator, is 
given by 



which leads to Eq. (J40 



The renormalization constant Z' 2 can also be obtained 
by computing Gn iQ assuming that the sources in Eq. 
(IA4[) are x- independent. In this case, by setting tq = 
we have 



G^( P )=p 2 + r 1 + (2u + el) 



+u 



-2e 2 



d d k 



(2?r) d k 2 + n 



d d k 



1 



2et 



d d k 



n -p 



(2n) d k 2 +r 2 J (2n) d [(p - k) 2 + r{\k 2 
^ 1 (A14) 



(27r) d k 2 ' 



n R = -±n. (A18) 
. z . . 

Since is finite, its dimension is just given by dimc- 
nional analysis, i.e., dim[n#] = d~ 2. Near the quantum 
critical point, we have that Z' 2 fZ z ~ /I 1 * 2 , where rj' 2 is the 
value of 72 [recall Eq. (|38|) ] at the infrared stable fixed 
point. Therefore, the scaling dimension of the Neel field 
is given by 



dim[n] = d — 2 — r( 2 - 



(A19) 



Thus, 



+2e 2 



1 - (2u + e 2 ) 
d d k 



d d k 



1 



(2ir) d {k 2 +n) 2 
1 



-2e 2 oP 2 



(27r) d [(p- k) 2 +n)k 2 
d d k 1 
(2^p[(p-fc) 2 + ri } 2 k 2 ' 



(A15) 



The anomalous dimension r/jv is defined by the scaling 
behavior (l32l). Thus, we also have that 



dim[n] 



and therefore we obtain 



d — 2 + r] N 



(A20) 



and 



G^ 2 (0,p;-p) = ^-G^ip) = ~u J ( 2 ?r ) d(fc2+ 1 r . 2)2 . 



Vn — d — 2 — 2rj' 2 . 



(A21) 



(A16) Using Eq. 



obtain the form given in Eq. ([5]). 
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